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ABSTRACT 

We combine new HST/ACS observations and existing data to investigate the wavelength depen- 
dence of near-IR (NIR) extinction. Previous studies suggest a power-law form for NIR extinction, 
with a "universal" value of the exponent, although some recent observations indicate that significant 
sight line-to-sight line variability may exist. We show that a power-law model for the NIR extinction 
provides an excellent fit to most extinction curves, but that the value of the power, /3, varies signifi- 
cantly from sight line-to-sight line. Therefore, it seems that a "universal NIR extinction law" is not 
possible. Instead, we find that as (3 decreases, R{V) = A{V)/ E{B — V) tends to increase, suggesting 
that NIR extinction curves which have been considered "peculiar" may, in fact, be typical for different 
R{V) values. 

We show that the power law parameters can depend the wavelength interval used to derive them, 
with the P increasing as longer wavelengths are included. This result implies that extrapolating power 
law fits to determine R{V) is unreliable. To avoid this problem, we adopt a different functional form 
for NIR extinction. This new form mimics a power law whose exponent increases with wavelength, 
has only 2 free parameters, can fit all of our curves over a longer wavelength baseline and to higher 
precision, and produces R(V) values which are consistent with independent estimates and commonly 
used methods for estimating R{V). Furthermore, unlike the power law model, it gives R{V)s that 
are independent of the wavelength interval used to derive them. It also suggests that the relation 
R{V) = -lM§^Evj ~ 0-79 can estimate R{V) to ±0.12. 

Finally, we use model extinction curves to show that our extinction curves are in accord with 
theoretical expectations, and demonstrate how large samples of observational quantities can provide 
useful constraints on the grain properties. 
Subject headings: ISMidust, extinction 



1. INTRODUCTION 

In the previous papers in this series, we began with 
an in-depth look at the structure and properties of 
UV extinction curves, inc l uding the 2 175 A "bump" 
(jFitzpatrick fc Massal Il986l . Il988l Il990f) . and most re- 
cently completed a general survey of Galactic extinc- 
tion co vering the near IR (NIR) th rough UV spectral 
regions (jFitzpatrick fc Massall2005al 2007, hereafter Pa- 
pers IV and V, respectively). These studies have uti- 
lized 2MASS JHK photometry, optical photometry, and 
low-resolution International Ultraviolet Explorer satellite 
{lUE) spectrophotometry. Our general goal of codifying 
the behavior of interstellar extinction over the widest 
possible wavelength range has been motivated by two 
scientific objectives: (1) to aid in the identification of 
the dust grain populations which produce the extinction 
and (2) to provide tools for the removal of the effects 
of wavelength-dependent extinction from astronomical 
data. 

During the course of this study, it has become clear 
that our ability to characterize extinction properties is 
strongly limited by the non-uniform type and quality 
of the data available in the various spectral domains. 
In fact, due to the success of the lUE mission, the 
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best characterized extinction curves are in the once- 
inaccessible UV region (1200-3000 A), while the worst- 
characterized regions are the "gaps" between UV and 
optical data ('^3000-3800 A) and between IR and opti- 
cal data ('^6000-10000 A). These regions are, in princi- 
ple, accessible from the ground but, due to calibration 
issues, remain poorly studied. In an attempt to rem- 
edy this situation, we obtained spectrophotometric data 
from the Advanced Camera for Surveys (ACS) aboard 
the Hubble Space Telescope (HST). The observations in- 
cluded both gap regions although, for reasons discussed 
in § 12.11 below, only the data spanning the optical-IR 
gap were useful for the program. These data, however, 
turn out to be of value not only for characterizing the 
wavelength-dependence of extinction along the chosen 
set of sight lines, but also for addressing the broader 
issue of a "universal" extinction law in the IR region. A 
variety of observations have suggested a common form 
to extinction at wa velengths from t he I band through 
^5 M m (e.g.. lRieke fc Lebofsky 1985: iMartin fc Whittetl 
Il990r ). consistent with a power-law form for the extinc- 
tion law, A\ (X , and an exponent value of P = 1.84 
for b oth diffuse and dark cloud lines of sight (jWhittet) 
l200 3f). More recent observations, however, have sug- 
gested significant sight line-to-sight line variability in IR 
extinctio n with a potentially large range of power-law ex- 
ponents (iLarson fc Whittetil2005H Nishivama et~aHl200d 
iFroebrich et al.ll2007l: iGosling et al. 2009). 

In this paper we use our new ACS observations in the 
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gap region between the optical and the NIR, along with 
existing data, to investigate the shape of the extinction 
law in this region and also the issue of whether a universal 
NIR law exists. In § 2 the target stars and the data used 
here (both new and archival) are described. In § 3, we 
discuss the technique used to derive extinction curves, 
and consider two families of analytic functions that can 
be used to represent them. In § 4, we examine how our 
results compare to theoretical models. § 5 summarizes 
the main results of the paper, and an appendix addresses 
technical issues concerning the ACS calibration. 

2. THE TARGET STARS AND THEIR DATA 

In this study, we examine extinction curves towards 
14 stars drawn from a set of more than 300 stars stud- 
ied in Paper V. These 14 stars, along with their spectral 
types, V magnitudes, and reddenings, are listed in Ta- 
ble 1. Since our method for deriving extinction curves 
involves modeling the observed spectral energy distribu- 
tions (SEDs) of reddened stars — using stellar atmo- 
sphere calculations and an analytical form of the extinc- 
tion curve — we require absolutely-calibrated photome- 
try and spectrophotometry. The observing program and 
the processing of the new HST/ ACS data are described 
in § 12. 1[ The ot her data sets utilized here (via the anal- 
ysis m IPaper Vl ) are described in § 12.21 and the complete 
energy distributions are presented in § 12.31 

2.1. New ACS Observations 

The new observations were part of our Cycle 14 
HST/ ACS program entitled "A SNAP Program to Ob- 
tain Complete Wavelength Coverage on Interstellar Ex- 
tinction." Our goal was to obtain ACS/HRC spectra in 
the near-UV (using the PR200L grating) and in the NIR 
(using the GR800L). These data, when combined with 
existing NIR {JHK) photometry, optical photometry, and 
lUE UV spectrophotometry, would provide complete ex- 
tinction curves over the range ~1150 A to ^2.2 /xm. We 
submitted a target list of 50 stars drawn from the sam- 
ple in Paper V. Ultimately, 22 stars were observed during 
2005 and 2006 and, of these, 14 were deemed suitable for 
this study. The ACS data were processed and calibrated 
using the aXe Spectral Extraction package (version 1.6), 
developed by M. Kiimmel, J. Walsh, and H. Kuntschner. 
We obtained the package from the Hubble Space Tele- 
scope European Homepage'^. The aXe software was run 
in conjunction with the IRAF/PYRAF package on a Sun 
Microsystems SunBladelOOO machine. 

Obtaining absolutely-calibrated spectrophotometry 
with the ACS is a challenging process in general and for 
our IIRC/PR200L observations it proved to be impossi- 
ble. The problem was the "red pile-up" effect whereby 
all photons between 4000 and 10000 A are focus ed onto 
a reg ion of the detector only 7 pixels wide (Larse iTet al.l 
l2006f) . This is particularly troublesome for red objects, 
such as our reddened targets, since diffraction spikes 
from the resultant bright spot overlap the blue region of 
the spectra. We found it impossible to extract reliably- 
uncontaminated spectra from our IIRC/PR200L obser- 
vations and ultimately excluded these data - and con- 
sideration of the UV-visual extinction gap ~ from this 
study. 

^ www.stecf.org / software / slitless^oftware /axe/ 



Processing the HRC/G800L observations was more 
straightforward and we followed the procedures outlined 
in the aXe Users Manual using an extraction slit height of 
±20 pixels. The result was a set of absolutely-calibrated 
first-order grism spectra for our 14 targets covering the 
range 5700-10200 A. Because absolute calibration is crit- 
ical to our program, we also verified the HRC/GSOOL cal- 
ibration by examining the HRC/GSOOL spectra of three 
DA white dwarf stars used to establish the calibration. 
We found a significant systematic discrepancy between 
the calibration targets processed by aXe and their corre- 
sponding model atmospheres. Consequently, we had to 
derive a correction factor for our spectra. This calibra- 
tion correction is provided in the Appendix. The final 
step in our ACS data processing was to trim the spectra 
to the wavelength range 6000-9500 A, where the fluxes 
are most reliable. 

The final set of 14 ACS/HRS/G800L spectra for our 
program stars are shown in Figure [TJ The spectra have 
been shifted vertically for clarit y. At the t op of the fig- 
ure, the SED of an ATLAS9 ( Kurucj|199l[) model atmo- 
sphere with Teff = 20000 K, log 5 = 4.0, and [m/H] = is 
shown for comparison. The model fluxes were computed 
m 20 A bins in the spectral range shown and smoothed 
by 2-points to approximate the 2-pixel resolution of the 
ACS data (i.e., ^^48 A). While our program stars span 
a relatively wide range in Toff (--15,000 - 40,000 K), the 
model SED shown provides a reasonable depiction of the 
stellar features expected in our stars. The H I Ha line 
at 6563 A and the upper H I Paschen lines near 9000 
A are the only spectral features of note for our targets 
and, at the resolution of the HRC/GSOOL, these features 
are very weak and generally undetectable. An exam- 
ination of the individual spectra show no signs of the 
H I lines although numerous low amplitude "bumps and 
wiggles" are present. These are instrumental features 
whose strength and location depend sensitively on the 
position of the stellar spectrum on the ACS detector. In 
the discussions to come, these features are smoothed out 
or eliminated so as not to influence the analysis. 

2.2. Other Data 

As will be described below in 3, we will use the 
SED-fitting results from iPaper VI . and thus implicitly 
rely on the data used there. These included low- 
resolution UV spectrophotometry from the lUE satel- 
lite, obtained via the Multimission Archive at STScI 
(MAST); ground-based [/BV^ photometry, obtained from 
the General Catalog of Photometric Data (GCPD) main- 
tained at the University of Geneva (Mermilliod, Mermil- 
liod, & Hauck 1997)^; and JHK photometry from the 
Two-Micron AU Sky Survey (2MASS), obtained from 
the 2MASS database at the NASA/IPAC Infrared Sci- 
ence Archive (IRSA)^. The only processing required for 
these data involved the lUE spectrophotometry. The 
lUE data at MAST were proces sed by the NEWSIPS 
software (Nichols & Linskv "1996"). As discussed in de- 
tail by Massa & Fitzpatrick ( 2000) , these data contain 
significant thermal and temporal dependencies and suf- 
fer from an incorrect absolute calibration. We corrected 

^ The GCPD catalog was accessed at 

|http://obswwwTun ige.ch / gcpd / gcpd . html 

° The ~2MASS data were accessed at 

|http://irsa.lpac.caltech.edu/applications/Gator| 
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the data for their systematic errors and pla ced them 
onto the HST/FOS flux scale of iBohlid (I1996D using the 
correc tion algorithms described by iMassa fc Fitzpatrickl 
(|2000f ). This step is essential for our program since our 
"comparison stars" used to derive extinction curves are 
stellar atmosphere models, so systematic errors in the 
absolute calibration would not cancel as they would for 
Pair Method of extinction determinations. 

2.3. A View of The Energy Distributions 

To provide a perspective on the SEDs used to derive 
the extinction curves in our studies, Figure [2] shows all of 
the SED data for our 14-star sample. The SEDs are ar- 
bitrarily shifted vertically for display. Small circles show 
the lUE spectrophotometry in the range A < 3000 A, 
and the ACS/HRC/G800L data in the range 6000 < 
A < 9500 A. 2MASS JHK values are shown as large 
filled circles in the NIK region (A > 10000 A). In the 
optical, UBV data are indicated as filled circles. To add 
to the coverage, we also show ground-based Stromgren 
uvby (filled triangles) and Geneva UB1B2V1G (filled di- 
amonds) photometry from the GCPD catalog whenever 
they are available. For all the photometry, the conversion 
to absolute fluxes was perform ed using the calibrations 
of iFitzpatrick fc Massal (|2005bl ). As noted above, most 
of the structure seen in the ACS/HRC/G800L data is 
instrumental. The dashed curves in the flgure are shown 
to guide the eye and to help isolate the individual stars. 
These lines are not model fits. 

3. THE ANALYSIS 

We begin this section by creating the extinction curves 
that will be used in the analysis. We then test whether 
the NIR extinction curves can be expressed as a simple 
power-law. We demonstrate that there is no "universal" 
exponent for the power law, and that the exponent which 
provides the best fit to a particular curve is sensitive to 
the wavelength range used to derive it. We also show that 
the R{V) determined from a power law fit to HD 164740 
is inconsistent with previous estimates and with an in- 
dependent estimate for R{V) based on its distance and 
absolute magnitude. These facts destroy our confidence 
in schemes which attempt to derive R(y) from power 
law fits, and motivate us to introduce an alternative, an- 
alytic form for NIR extinction curves. We show that this 
new form actually describes the NIR data as well as any 
power law, and that it does so over a much larger wave- 
length range. Further, the new formulation produces an 
R{y) for HD 164740 which is consistent with the inde- 
pendent estimate. It also gives R{V) values for the other 
program stars that agree with previous estimates. 

3.1. THE NIR EXTINCTION CURVES 

We produced normalized NIR extinction curves, k{\ — 
V) = E[\-V)/E{B-V) for our 14 sight lines by adopt- 
ing the a ppropriate stellar atmosphere models de rived in 
iPaper Vl as the unreddened SEDs. In iPaper Vl the ob- 
served SEDs of 328 reddened stars — including those 
studied here — were fit with a combination of stellar at- 
mospheres (to estimate the unreddened SED) and a flexi- 
ble analytical representation of interstellar extinction (to 
model fc(A — y), utilizing a uniform set of lUE UV spec- 
trophotometry, UBV optical photometry, and 2MASS 
JHK NIR photometry. 



For each progr am star h ere, we adopt the best-fit in- 
trinsic SED from lPaper Vl . which is defined by the values 
of Toff, logg, [m/H], and Vturb listed in their Table 3. The 
extinction in the JHK bands was computed using syn- 
th etic photo metry of the stellar model fluxes as described 
in lPaper Vl and references therein. To mitigate noise and 
more clearly view trends in the ASC data, we rebinned 
the ACS measurements and model SEDs into 500 A bins 
centered at 6250, 6750, 7250, 7750, 8250, 8750, and 
9250 A, thereby utilizing all of the ACS data within the 
6000-9500 A range of our measurements. These fluxes 
were then converted into 7 independent magnitudes. To 
smooth instrumental features in the ACS data, we fit- 
ted the individual spectra with fourth-order polynomials, 
employing a sigma-clipping algorithm. The results were 
smooth representations of the ACS data, from which the 
synthetic magnitudes were computed. We list the values 
of the normalized NIR extinction for the program stars 
in Table 2. The labels ''ACS62;' ''ACS67;' etc. refer to 
the seven synthetic filters described above. 

The l-cr errors listed in the table incorporate both 
the uncertainty in the intrinsic SEDs and the uncertain- 
ties in the ACS, JHK, and F-band measurements from 
which k{X — V) is computed. For each star, we used the 
100 Monte Carlo error simulations described in Paper 
V, which each yielded an independent estimate of the 
shape of the intrinsic SED. For each of these SED's, we 
computed 20 NIR extinction curves, each using the ob- 
served photometry convolved with a random realization 
of the expected uncertainties in the measurements. Our 
noise model includes a 1-a error of ±0.015 mag in V and 
the standard errors provided with the 2MASS JHK pho- 
tometry. However, while calibrati ng synthetic photom- 
etry fo r use in SED modeling, we (jFitzpatrick fc Massal 
I2005bf) found that the 2MASS magnitudes - notably H - 
had a larger scatter than could be explained by the stated 
uncertainties. We thus quadratically combine the stan- 
dard 2MASS uncertainties with values of 0.007, 0.040, 
and 0.017 mag for J, H, and K, respectively, to arrive 
at the total uncertainty for each magnitude. Finally, we 
adopted the point-to-point photometric errors produced 
by the aXe processing software for the ACS/HRC/G800L 
spectra and a 1-a zero-point uncertainty of ±2.5% (see 
the Appendix). The final result was an ensemble of 2000 
curves for each star, representing the likely range of un- 
certainties. The errors listed in Table 2 are simply the 
standard deviations of each point derived from the 2000- 
curve set. 

The data from Table 2 are plotted in Figure[3l with the 
key to the symbols given in the legend. Three curves are 
over plotted for most of the stars. These are power-law 
model fits which will be discussed below. 

3.2. THE POWER-LAW MODEL 

We begin by comparing our extinction curves to a basic 
power-law extinction model of the form 

^(A) cx A^'' , (1) 

where ^(A) is the total extinction at wavelength A and (3 
is a wavelength-independent exponent. Converting this 
to the normalized form of our extinction curves yields 

MA -v)^ e{b'_v] = ~ ^^^^ ' 
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where k{X — V) is a function of the power /3, a scale 
factor kiR, and the well-known ratio of selective to total 
extinction RiV) [= A{V)/E{B - V)]. In the following 
discussion, we use non-linear minimization routines 
(jMarkwardtl 12009) to fit the observed extinction curves 
with equation ^ for two different cases. In the first, km 
and R{V) are taken as free parameters and a "universal" 
value of /3 — 1.84 is assumed (the "/3-Fixed" case). In 
the second, km, (3, and R{V) are all allowed to vary 
(the Variable" case). Previous studies have shown 
that when extinction curves are modeled by equation ([1]) , 
the fits tend to diverge from observations at wavelengths 
shortward of the / photometric ban d, i.e., at A < 7500 
A (se e, e.g., Figures 1, 2, and 3 of iMartin fc Whittetl 
|1990[) . Likewise, we found that, for most of our sight 
lines, the inclusion of ACS data shortward of 7500 A 
greatly degraded the quality of the fits by increasing x^- 
Thus, we only use data with A > 7500 A in our power 
law analyses. These data are shown in Figure [3] as large 
filled symbols. 

The best-fits to our data using equation ^ are shown 
in Figure [3] as the dashed ( "/3-Fixed" ) and solid ( "/3- 
Variable" ) curves in each panel. The curves are extrap- 
olated over the full wavelength extent of the panels, but 
were determined using only the data indicated by the 
large filled symbols. In most of the panels, it is very 
clear that the power-law functions fail badly for A < 7500 
A — which led to the wavelength restriction described 
above. The dotted curves shown in the Figures are /3- 
Variable fits to the ACS data alone, utilizing all seven of 
the synthetic ACS magnitudes. These curves address the 
ability of a specific parameterization to provide a reason- 
able and consistent extrapolation to longer wavelengths. 
This important issue will be discussed further below. 

Table 3 presents the numerical results of the fits illus- 
trated in Figure 3. The 1-a uncertainties listed in the 
table are determined from the ensemble of 2000 Monte 
Carlo noise model realizations computed for each sight 
line. For each parameter they represent the dispersion of 
the results from power-law fits to the simulated curves 
about the best-fit results. The final column of Table 3 
gives the F-statistic, constructed from the x^ values for 
the /3-Fixed and /3-Variable cases. The F-statistic can 
be used to determine whether the addition of a new fit- 
ting term (in this case /3) significantly improves the fit. 
Large F values, i.e., F > 10, indicate the vahdity of the 
new parameter while progressively smaller F values indi- 
cate progressively less justification (see, e.g.. iBevingtonI 
[1969). Although, to the eye, most of the fits in Figure 3 
look reasonable, the F-statistic demonstrates that the 
fits for at least half of the sight lines are significantly im- 
proved by allowing a variable /3, strongly justifying the 
inclusion of /3 as a free parameter. The largest F values 
tend to occur for the sight lines whose R{V) values differ 
most from the Galactic mean of 3.1, particularly those 
with the largest RiV). The inability of the /3 = 1.84 
model to adequately represent large R{V) sight lines is 
readily apparent in the last panel of Figure [H 

An essentially identical method of assessing the sta- 
tistical significance of the derived /3 values is simply to 
compare them with the "universal" value of 1.84. If the 
new values differ from 1.84 by 2-cr or more — which is 
the case for seven of our sight lines — then there is very 



strong statistical justification for rejecting 1.84 as the 
appropriate value of /3. A visual illustration of the sig- 
nificance of the derived /3's is shown in the top panel 
of Figure IH where we plot /3 from the /3- Variable case 
against E{B — V) for each sight line. The dotted line is 
/3 = 1.84. The sample mean and median values of /3 are 
1.78 and 2.00, respectively. These are not far from the 
universal value, although the figure clearly illustrates sig- 
nificant departures exist. The figure also demonstrates 
that 1.84 cannot be replaced by some better-determined 
universal value. Indeed, it challenges the notion of the 
existence of a universal IR extinction law. 

A further challenge to the universal extinction model 
is shown in the lower panel of Figure |31 where we plot /3 
against R{V) for each sight line. Note that these R{V) 
values are determined by extrapolating the power-law 
models to A^^ = — a practice which will be dis- 
cussed further below. As might be anticipated from equa- 
tion ([2]), the uncertainties in /3 and R{V) are correlated, 
and the directions and strengths of the correlations are 
indicated in the figure by the l-a error bars. For com- 
parison, we als o show the values o f /3 an d R{V) for the 
star p Oph from lMartin fc Whittetl (|1990l filled triangle). 
The consistency of the power-law exponent for this dense 
cloud sight line with the results for diffuse cloud sight 
lines has been a significant motivation for the percep- 
tion of a universal P ~ 1.8 law (e.g.. Fig. 3.7 in IWhittetl 
|2003( ). Although the p Oph measurements were not made 
in the same manner as ours, they utilize data over the 
same wavelength range (/ through K-hand) and should 
be compatible. The data in Figure [3] reveal a striking sys- 
tematic trend — /3 is anti-correlated with R{V) in the 
sense that la,ige-R{V) sight lines have smaller than aver- 
age values of (3. As R{V) is generally considered a coarse 
indicator of grain size, i.e., populations of larger-than- 
avcragc grain s yield larger-t han-average R{V) (e.g., see 
iMathis fc WallenhorstJ [l98l[ ). our results indicate that 
the shape of IR extinction curves parameterized by /3 
are environmentally-dependent. As shown in Paper V, 
there is considerable intrinsic scatter in correlations be- 
tween R{V) and other properties of interstellar extinc- 
tion curves, and often the correlations are only apparent 
when the most extreme R{V) values are considered. Our 
results suggest that previous analyses have not had ac- 
cess to a wide-enough range in R{V) to reveal the envi- 
ronmental dependence of /3. 

Since a major motivation for using an analytic form 
for extinction is to extrapolate the curve to infinite wave- 
length and estimate R{V), it is troubling that the R{V)s 
derived from the power law curves depend strongly on the 
adopted exponent. Figure 3 and Table 3 show that the 
/?-Fixed and /3- Variable approaches give similar results 
when R{V) is near the Galactic mean of '^3.1, but yield 
increasingly divergent results at both larger and smaller 
R{V). This is best illustrated by the HD 164740 sight 
line. This is the most extreme sight line in our sample 
in the sense that it has the smallest /3 and largest R(V). 
Furthermore, there is no indication that it is pathologi- 
cal, since Figure H] shows that it is consistent with an ex- 
trapolation of the trends exhibited by less extreme sight 
lines. The problem with HD 164740 is that, even though 
the /3- Variable fit provides a fine representation for the / 
through if -band extinction (see Figure [3|), the extrapo- 
lation of its power-law fit to longer wavelengths is incon- 
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sistent with the physical expectation that '^^^^^ 







as A ^ cxD, (which demands /3 > 1). 

An additional, but related, problem is that the R{V) 
determined from the /3-Variable fit for HD 164740 is 
demonstrably incorrect. HD 164740 is a member of the 
cluster NGC 6530, whose distance, d, has been deter- 
mined independently of HD 164740. Therefore, given its 
absolute magnitude. My, it is possible to derive an in- 
dependent estimate for R{V) by rearranging the formula 
for the distance modulus. Specifically, 

y-Afv-51og(d/10pc) 
= E{B-V) ' 

The data needed to compute R{V) for HD 164740 are 
d = 1330 pc (Paper V), My = -4.8 (lWalbornlll973h . V 
and E{B — V) (see Table 1). We consider two cases which 
span the range of expected absolute magnitudes: a single 
star or a binary with identical components. Using these 
data for HD 164740, we find 5.2 < R{V) < 6.1. These 
values are significantly smaller than the R{V) — 7.0 in- 
ferred from extrapolation of the /?- Variable power-law fit. 
One can appeal to uncertainties in the distance and abso- 
lute magnitude of HD 164740, but these are insufficient 
to explain such a large discre pancy. In support of this 
conclusion. iHecht et al.l (|1982f) find R{V) = 5.6 based on 
a comparison of the brightness and color of HD 164740 
with those of the nearby - and presumably equally dis- 
tant — 0-type star HD 164794 (9 Sgr). Note that other 
estimates of R{V) for HD 164740 have been made by ex- 
trapolating an assumed form for NIR extinction by, for 
example, Cardclli et al. ()1989l ) [R{V) = 5.30] and more 
recently by us in Paper V [R{V) — 5.21]. However, the 
assumed forms, which are (or resemble) power-laws with 
large exponents {/3 =1.6-1.84) are clearly inconsistent 
with the observed shape of the NIR extinction curve of 
HD 164740. Thus, these results cannot be used to prove 
an error in the /3- Variable value. 

The HD 164740 discussion shows that, even though 
a power-law can fit an NIR extinction curve quite well 
over a specific wavelength interval, extrapolation of such 
fits to longer wavelengths can be unre liable. This was 
anticipated bv lMartin fc Whittetl (|1990f ) who noted that 
it would be "remarkable" for a single power-law to be 
appropriate over an extended wavelength range. Never- 
theless, this caveat is explicitly ignored when a single- 
exponent power-law is used to estimate R(V). To em- 
phasize this point, consider the dotted curves in Figure 3. 
These are power- law fits (eq. [2]) to all of the ACS data. 
I.e., covering the full 6000 - 9500 A wavelength range. 
While they certainly provide an excellent analytical rep- 
resentation of fc(A — V) at these wavelengths (x^ < 0.2 
and, in 12-of-14 cases, < 0.1), nevertheless it is quite 
clear that extrapolations of these fits to longer wave- 
lengths are grossly inconsistent with the JHK data and, 
consequently, the intercepts of the extrapolated curves 
are totally unrelated to R{V). In addition, the /3s are 
always smaller than those inferred from the / through 
-ftT-band fits in Table 3 and are nearly always less than 
the limiting value of one. These results suggest a system- 
atic effect in which the exponent of a best-fit power-law 
systematically decreases as shorter wavelengths are con- 
sidered and the implied R{V) increases. 

At this point, several straightforward conclusions can 
be drawn. 



1. All of our NIR extinction curves can indeed be well- 
represented by a simple power-law formulation in 
the J-band through K band spectral region. 

2. However, the power-law exponent, /?, varies signif- 
icantly from sight line-to-sight line, and is incom- 
patible with the thesis of a universal form to NIR 
extinction. 

3. The value of the exponent, /3, is anti-correlated 
with R{V), suggesting that the detailed shape of 
NIR extinction curves is likely a function of the 
line of sight grain size distribution. 

4. The parameters of power law fits can depend quite 
strongly on the wavelength interval used to deter- 
mine them, with the value of /3 tending to decrease 
as shorter wavelengths are included. This result 
implies that extrapolations of power law fits to 
longer wavelengths can be an unreliable means for 
estimating R{V). 

3.3. AN ALTERNATIVE MODEL 

Given the problems encountered with the simple power 
law representation for NIR extinction, we decided to 
search for a different functional form, whose parame- 
ters can be estimated from accessible data and can be 
extrapolated to long wavelengths to obtain consistent, 
meaningful estimates of R{V). W hile many such forms 
are possible, the one introduced bv lPeil (|1992l ) in a study 
of Milky Way and Magellanic Cloud extinction produced 
remarkably good results^. We use a generalized version 
of the function, with more degrees of freedom. It has the 
form 



k{X-V) = k 



1 



IR. 



1 + (A/Ao 



^ - RiV) , (4) 



where a, /c/^, Aq and R{V) are potentially free param- 
eters. For A ^ Aq, this function reduces to a power law 
with exponent —a. For A Ao, it flattens to a constant 
whose value is km — R{V). Between these extremes, 
it resembles power law whose exponent increases with 
wavelength — which is exactly the behavior suggested 
by our analysis in the previous section. 

While the introduction of four free parameters may 
seem a steep price to pay for abandoning the simpler 
power-law formulation, we quickly discovered that two 
strong observational constraints can be placed on the pa- 
rameters. The first is that Aq can be replaced by constant 
for all of our sight lines, which span a wide range of ex- 
tinction properties, without significant loss of accuracy. 
The second is that the scale factor, kj^, can be expressed 
as a linear function of R{V), again without loss of accu- 
racy. These constraints result in a simplified version of 
equation ^ given by 



kiX-V) = [0.349-1-2. 087i?(y)] 



1 



1 + (A/0.507)" 



RiV) , 
(5) 



^ We thank Karl Gordon for bringing this form to our atten- 
tion. It was used in a study of Galactic extinction curves by 
IClavton et"all (200% an d a m ore general discussion of its utility 
will be given in IGordonI | |2009|) 



6 



Fitzpatrick & Massa 



where A is in fim. Equation (O contains only two free 
parameters, a and R{V); the same number as the f3- 
Fixed power law model, and one less than the /3- Variable 
model. 

Equation ^ proved so successful in reproducing the 
shapes of the NIR extinction curves that we applied it 
to the JHK data and the entire wavelength range of the 
ACS data. The results are given in Table 4 and Figures 5 
and[ni In Figure 5, the solid portion of the curves indicate 
the wavelength interval over which the fit was performed, 
and the dotted portions show extrapolations of the fit to 
longer and shorter wavelengths. Note that, although the 
y-band point was not included in the fits, extrapolations 
of them pass quite close to it. In fact, the V-band could 
have been included in the fits without much loss of overall 
accuracy. Table 4 gives the a and R{V) that result from 
the fits and their uncertainties. These are plotted in 
Figure [6l The correlated uncertainties were computed as 
described in § 13.21 The range in a and the relationship 
between a and R{V) mimic those of the /3- Variable fits 
discussed in § 13.21 

That equation ([5]) provides a much better fit to the ob- 
served NIR extinction curves can be seen by comparing 
the mean values listed in Tables 3 and 4 for the dif- 
ferent models. The mean values for the /3-Fixed, and 
equation ([5]) are 0.78 and 0.24, respectively. Thus, even 
though the fits using equation ([S]) include more points, 
their mean is much better than the 2-parameter f3- 
Fixed power law model. Furthermore, unlike the /9-Fixed 
model, equation ([5]) provides a good fit for all of the 
curves. 

Notice that although the mean x^ for the 3-parameter 
/?- Variable from Table 3 is slightly better than the equa- 
tion ([5]) value (0.17 compared to 0.24), a fairer compari- 
son between the two is given in Table 4. The last column 
of the table lists the values that result from fitting f3- 
Variable models to the same data as the equation ([5]) fits. 
In this case, we see that the mean for the /3- Variable 
fits is 0.44, considerably larger than the 0.24 that results 
from the equation ([5]) fits. 

An additional advantage to the formulation in Equa- 
tion ([5]) is its robustness to changes in the wavelength 
coverage of the data. For instance, if we fit only the ACS 
data, spanning the wavelength range 6000 - 9500 A, then 
the resultant values of R{V) differ on average by only 
-1-0.09 from the results in Table 4, with a mean scatter 
of ±0.27 and no systematic trends. Likewise the implied 
a's differ on average from those in Table 4 by —0.05, 
with a scatter of ±0.20. Thus, in striking contrast to 
the /3- Variable reliable and stable estimate of IR 

extinction can be derived, even in the absence of longer 
wavelength JHK data. 

It is important to determine whether the values of 
R{V) determined by equation ^ agree with previous 
estimates. To begin, we note that the R{V) value deter- 
mined for HD 164740 is 6.00, within the expected range 
described in § 13.21 We can also make a broader compar- 
ison by comparing the values in Table 4 to a commonly 
used means to estimate R{V), which is based on optical 
and JHK photometry. This approach estimates R{V) 

from the relation R(y) = —'^■^ e[b-v) ' which is based 
on v an de Hulst's theoretical extinction curve No. 15 
(e.g.. lJohnsonlfl968l ). While, in general, our results show 



that two parameters are needed to fully define an NIR 
extinction curve, this does not exclude the possibility 
that relationships may exist between specific properties 
of the curves, at least to a high degree of approximation. 

Figure [5] is a plot of ^I^Zyj (Table 2) versus R{V) from 
Table 4. The two quantities are clearly strongly corre- 
lated. The van de Hulst-based formula (dashed line) 
follows the general trend for low and moderate values 
of R{V), but systematically underestimates it for large- 
R{V) sight lines. Nevertheless, the overall agreement is 
quite good, with an RMS deviation of 0.26 in R{V), im- 
plying that the R{V) values determined by equation (O 
are in general agreement with previous estimates. 

In fact, if we assume that the R{V) derived from equa- 
tion ([5]) are correct, we can derive a more accurate means 
for determining R{V) from the broad band photometry. 
The solid line in Figure [8] is the unweighted linear fit, 

R(V) = ~1-^G e[bZv) - 0-79, which has an RMS devia- 
tion for R{V) of 0.13. Note that both lines yield similar 
resuhs for R{V) < 3.5. 
The major conclusions from this section are: 

1. Equation ^ mimics a power law whose exponent 
increases with wavelength, which is exactly the be- 
havior our analysis in 13.21 implied was needed. 

2. Initial fits to the data demonstrated that the num- 
ber of free parameters needed to fit NIR curves can 
be reduced to 2 (see, eq. [5]). 

3. Equation ([5]) provides good fits to all of the ex- 
tinction curves, and even performs better than the 
/3- Variable model, which has one more free param- 
eter. 

4. Equation (O produces R{V) values which are con- 
sistent with independent estimates, and are in ac- 
cord with commonly used methods for estimating 
R{V). 

5. Equation ([5]) gives R{V) values that do not depend 
strongly or systematically on the wavelength inter- 
val used to derive them. 

6. For our data, the relation R{V) = — 1.36 ^|^~yj — 

0.79 reproduces R(V) to within ±0.12 over a large 
range in R{V). 

4. DISCUSSION 

In general, there are two motivations for determining 
an analytical model for the shapes of interstellar extinc- 
tion curves. The first is that it allows a consistent in- 
terpolation or extrapolation of extinction properties into 
spectral regions where measurements are lacking. This 
is used to define R{V) and, hence, A{V). The second is 
that the functional dependence of extinction over some 
particular spectral region may shed light on the physical 
processes controlling the extinction, e.g., the composition 
and size distribution of the grains responsible for the ex- 
tinction. Further, the analytic representation may yield 
specific parameters whose range of variability could be 
related to how the grain populations respond to different 
physical environments along the line of sight. The first of 
these goals was addressed in the previous section. Here, 
we focus on the second. 
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Our results suggest that at least two parameters are 
needed to characterize the NIR extinction, R{V) and a. 
It seems reasonable to suggest that R{V) provides a mea- 
sure of the mean grain size along the line of sight, and 
that a is probably related to the grain size distribution. 
This is because, in the small particle limit (where the ra- 
tio of grain-size to wavelength is small - as in the IR), Mie 
scattering theory implies that both absorption and scat- 
tering are expected to approach power law forms with 
an exponent of 4 for sc attering and ^^1 for absorption 
(see Chapter 7 of ISpitz cr 1978). Power-law exponents 
determined in the NIR generally lie between these two 
extremes. Given the sensitivity of extinction to grain 
properties, however, it would be surprising for extinction 
to adhere strictly to a single-exponent power law form 
over an extended wavelength range (Martin fc Whittetl 
fl990l) . 

In fact, these arguments make the suggestion of a uni- 
versal power-law exponent for the IR troubling. It is 
clear from the study of other aspects of extinction (e.g., 
the shape of UV extinction curves) that the properties of 
grain populations differ significantly from p l ace to place 
in interstellar space. As iMartin fc Whittetl (iQQfih note, 
since the absorption and scattering properties of grains 
are generally strongly dependent on the grain compo- 
sition and size, the suggestion of universality requires 
highly constrained sight line-to-sight line modifications 
in grain properties which must occur in such a way as 
to yield, for example, a range in R{V) but a single IR 
power- law exponent. Our results remove this troubling 
constraint. We find that the exponent is spatially vari- 
able and likely does reflect underlying differences in the 
grain populations. 

To verify that the observational parameters obtained 
from equation (O are consistent with grain models and 
to determine whether they can be used to constrain 
the grain parameters, we compa re our results with the 
theore tical calculations given by iMathis fc WallenhorstI 
()1981l hereafter MW81). Alt hough m ore recent and so- 
phisticated models exist, the iMWSlI models are simple, 
accessible, and adequate to demon strate th e points we 
wish to make. In their Table 1, iMWSlI present the 
extinction optical depths (per 10^^ H atoms) for sep- 
arate populations of silicate and graphite grains, each 
wi th a power-law siz e distribution (following the work 
of IMathis et al.lll977l hereafter MRN) and six different 
combinations of the upper and lower grain-size cutoffs, 
a+ and a_, respectively. We produced model extinction 
curves by summing the graphite and silicate extinction 
for all possible combinations of the six size distributions 
(of which there are 36) and three different assumed val- 
ues for the dust-phase C/Si ratio. As a baseline, we 
adopt a ratio {C/Si)dust = 6.5, which is consistent with 
the mean for the variety of dust grain models considered 
by IZubko et al.l (|2004l ) and also with mean dust-phase 
abun dances based o n solar composition (see, e.g., Table 
2.2 of lWhittet]l2003[ ). We used this value for one set of 
36 model optical depth curves and values a factor-of-2 
lower (3.25) and higher (13) to produce two additional 
sets. After constructing this ensemble of 108 curves, we 
converted them from optical depth to absolute extinc- 
tion and then normalized them by E{B — V), producing 
curves of A{X)/E{B — V). The justification for producing 
this multitude of curves is simple: we have no a priori 



reason to believe that the grain-size distributions must be 
identical for the graphite and silicate grain populations, 
as the destruction and, certainly, the growth mechanisms 
are unlikely to be identical. Furthermore, we do not know 
the exact abundances of C and Si in the ISM, and suspect 
that the dust-phase C/Si ratio varies from sight line-to- 
sight line. Thus our model curves represent a wide range 
of conditions that may occur in interstellar space. 

We seek to compare the shapes of model IR curves 
with specific R{V) values to the observed parameters de- 
termined from equation ([5]). Since the iMW81i curves are 
presented at only at a small number of wavelength points, 
including optical/IR points at B, V, 1.11 fim~^ and 0.29 
/im~\ we take R{V) as the value of A{X)/E{B - V) at 
the V point and characterize the shape of the IR with a 
power law fitted through the points at 1.11 /im~^ (A — 
0.9 /zm) and 0.29 fj,mr^ (A = 3.4 /im), yielding an expo- 
nent we refer to as "aMWsi-" This provides a measure 
of the general curvature of the extinction between and 
1.11 /im^ -1 (i.e., A > 9000A). A fair comparison with the 
predictions of equation ^ requires that we compute a 
similar value oiaMWSi for each of the curves in Figure 5, 
using equation ([5]) at the same wavelengths. 

The comparison between aMWsi and R{V) for the 
models and the observations is shown in Figure [T] Filled 
circles show the observational results. The distribution 
of points is similar to that in Fig. [6l except that the 
values of aMWSi are ~ 0.15-0.20 smaller than a. Mea- 
surement errors were determined in the same way as de- 
scribed above. The theoretical aMWSi and R{V) values 
derived from the 108 model curves are depicted as the 
x's. The symbol sizes correspond to the (C/Si) ratio in 
the dust, with {C/Si)dust — 3.25, 6.5 and 13 represented 
by small, medium, and large crosses, respectively. Two 
interesting points are evident in the Figure. First, the 
main distribution of model points at large R{V) aligns 
nicely with the "observed" values of R{V) and aMwsi- 
This indicates that the functional form proposed in equa- 
tion (O is compatible with theoretical calculations. Sec- 
ond, there is a group of model points that do not over- 
lap the domain covered by the observations (the models 
which have small aMwsi o,'<^d small R{V) values). This 
may indicate that there are specific size distribution and 
abundance combinations which are not favored by na- 
ture. In particular, most of these points arise from the 
combination of a silicate grain population with a small 
upper size cutoff (a+ — 0.25) and a graphite grain popu- 
lation with a large upper si ze cutoff (a+ = 0.40 or 0.50). 
Note that, because IMWSII were explicitly interested in 
studying the behavior of curves with large R{V) values, 
they considere d grain size distributions with larger mean 
sizes than the 'MRN distribution, which provides a rea- 
sonable fit for the Galactic mean curve with R{V) = 3.1. 
As a result, there are no models with large (3 and small 
R{V) to compare to the observations. Additional cal- 
culations and mor e sophisticat e d gra in models, such as 
those discussed bv' Zubko et al.l ([2004) , will be needed to 
address this. However, it is the large R{V) sight hues 
where our results differ most strongly from the paradigm 
of a universal power law and, in these cases, the compar- 
ison in Figure [7] is quite satisfactory. 

From the preceding discussion, we conclude that the 
two parameters needed to fit the observations using equa- 
tion 15]) are probably the minimum required. This means 
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that the hope of finding a one parameter representation, 
such as a universal power law with an adjustable R{V) o r 
a model of the sort described by ICardelli et al.1 ()1989D . 
where the curve shape can also have an R{V) depen- 
dence, are probably unjustified. We were also able to 
demonstrate how examining the observed parameters can 
provide useful constraints on the physical parameters of 
dust models. 

5. SUMMARY AND FINAL COMMENTS 

Collecting the results from the previous sections, we 
have drawn the following conclusions from our analysis 
of NIR interstellar extinction along 14 sight lines: 

1. Our NIR extinction curves can be well-represented 
by a power-law formulation in the /-band through 
K band spectral region. However, the exponent of 
the power-law varies significantly from sight line- 
to-sight line, which is incompatible with the notion 
of a universal form to NIR extinction. 

2. The value of exponent in the power law, /?, is anti- 
correlated with R{V), showing that the shape of 
NIR extinction curves are likely functions of the 
line of sight grain size distributions. 

3. The parameters derived from a power law fit can 
depend quite strongly on the wavelength interval 
used to determine them, with the value of (3 in- 
ferred from the fit increasing as longer wavelengths 
are included. This result demonstrates that extrap- 
olations of power law fits to longer wavelengths are 
an unreliable means for estimating R{V). 

4. Because of the biases in the power law, we have 
adopted a different form for the NIR extinction, 
equation ([4]), which has the following properties: 

• Equation ([4]) mimics a power law whose expo- 
nent increases with wavelength - a property 
demanded by the power law analysis. 

• The 4 free parameters in equation ([4]) can be 
reduced to 2 (see, eq. [5]), without affecting 
the quality of the fits to the NIR extinction 
curves. 

• For R{V) ~ 3, the curve produced by equa- 
tion ([5]) is very similar to the commonly used 
/3 = 1.84 power law. 

• Equation ([5]), with its 2 free parameters, pro- 
vides better fits to all of the extinction curves 
than the power law model which allows the 
exponent to be a free parameter - giving it 3 
free parameters. 

• Equation (O produces R{V) values which are 
consistent with independent estimates, and 
are in accord with commonly used methods 
for estimating R{V). 

• Unlike the power law models, equation ([5]) 
gives R{V) values that do not depend strongly 
or systematically on the wavelength interval 
used to derive them. 



5. For our data, the simple relation R{V) ~ 
— 1.36 ^1^2^1 — 0.79 reproduces the R{V) deter- 
mined from equation ([S]) to within ±0.12, over a 
large range in R{V). This is a useful relation for es- 
timating R{V) when only limited broad band pho- 
tometry is available. 

6. Using theoretical grain models, we showed that the 
two parameters needed to fit the observed curves 
using equation ([5|) are probably the minimum re- 
quired. As a consequence, the notion of finding 
a one parameter representation, such as a uni- 
versal power l aw w ith an adjustable R{V) or a 
ICardelli et al.l ()1989) -like model where the curve 
shape can also have an R{V) dependence, are prob- 
ably unjustified. 

7. We demonstrated how examining the observed pa- 
rameters can provide useful constraints on the 
physical parameters of dust models. 

Our results bear heavily on the recent obser- 
vatio ns of "n on-univer saF' /3 values n oted in § 1 
Ci.e.. iLarson fc Whittc t] l2005t iN ishiva mTerall 120061 : 
iFroebrich et~aITl2007i: iGosling et all i2009f ) and indicate 
that they should be viewed with some caution. On the 
one hand, these may reflect the intrinsic variability ex- 
pected from an interstellar medium with a wide range in 
physical properties and grain processing histories for the 
specific conditions along the respective sight lines. How- 
ever, care must be taken in comparing [3 values from dif- 
ferent datasets. Our results show that power-law models 
for a given sight line will yield different exponents, de- 
pending on the wavelength range used to derive them. 
For example, a Variable power-law fit to equation [5] 
in the JHK region gives exponents that are ^^0.3 larger 
than those determined by fits which also include the I- 
band. This effect is likely present to some degree in the 
above studies, which are based solely on JHK data, and 
should be taken into account before we can obtain an 
accurate picture of the range in IR extinction properties. 

Finally, we point out that while our results are inter- 
esting, they are based on only 14 sight lines. Clearly, 
more NIR data are needed to verify the generality of the 
constraints placed on equation (jj]) to obtain equation ([5]), 

and to verify the relationship R{V) = — 1.36 ^|g~^j — 
0.79. 
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APPENDIX 

ABSOLUTE CALIBRATION OF THE ACS/HRC/G800L DATA 

Because the absolute calibration is critical to our program, we took the additional step to verify the HRC/G800L 
calibration by processing G800L spectra for three spectrophotometric standard DA w hite dwarfs (G 191B2B, GD 71, 
and GD 153) used to help establish the fundamental HST Rux cahbration (see, e.g.. lBohIiiilll996D . Our procedure 
was simple. We downloaded ACS/HRC/G800L spectra for these stars from MAST, processed them in exactly the 
same manner as our target stars, and then compared the resultant spectra with the model fluxes used to establish 
the absolute calibration. The models were obtained from the CALSPEC Calibration Database at the Space Telescope 
Science Institute"^ (the files gl91b2bjnod_004.fits, gd71.mod.005.fits, and gdl53jnod.004.fits). 

Figure lAll shows the results of this comparison. The different symbols show the ratio of the model fluxes to the 
ACS fluxes as a function of wavelength for each of the three stars. Clearly, the processed IIRC/G800L fluxes are 
increasingly underestimated at wavelengths longward of ~6000 A, with the effect increasing dramatically below ^9500 
A. Fortunately, the flux discrepancy is systematic and consistent among the three standards and, thus, can be corrected 
in a simple way. We derived a correction curve for our processed spectra by fitting a polynomial to the data in Figure lAT] 
and the result is shown as the smooth solid curve in the figure. The systematic scatter of the three stars about the mean 
correction suggests a zero point stability in the fluxes of ~2.5% at wavelengths shortward of 9500 A. The variability 
increases significantly at longer wavelengths. The mean correction in Figure lAll is applied as a multiplicative factor 
to the calibrated results of the aXe processing. 

To perform an independent test of our calibration correction, we examined its effect on the IIRC/G800L spectrum 
of the Galactic 07 star HD 47839 (15 Mon). ACS spectra for HD 47839 were obtained from MAST, and processed in 
the same way as our program and spectrophotometric calibration stars. Figure [X2] shows the multi- wavelength SED of 
HD 47839 spanning the range -1150 A to —2.2 /im. In the main figure are lUE low-resolution UV spectrophotometry 
(A < 3000 A), Johnson, Stromgren, and Geneva optical photometry (3300 < A < 6500 A), and Johnson i?/and 2MASS 
JHK NIR photometry (A > 6500 A). The sofid curve is a 37,000 K TLUSTY s tellar atmosphere model (iLanz fc Hubenvl 
l2003f ) fitted to the data in the manner described in our previous papers (^e.g.. lFitzpatrick fc Massall2005al lbl). The inset 
to the figure illustrates the effect of our HRC/G800L correction curve. The smooth curve shows the best fit model 
while the corrected and uncorrected G800L data are over plotted. The HD 47839 data clearly confirm the results from 
the analysis of the spectrophotometric standards and the effectiveness of the calibration correction in Figure lAll At 
wavelengths shortward of 9500 A the wavelength-dependent signature of the calibration deficiency has been removed 
and the mean flux level lies about 1% above that predicted by the best-fit SED model, consistent with the scatter 
among the spectrophotometric standards. At wavelengths longward of 9500 A, the fluxes are clearly much less reliable. 

The calibration correction curve derived in Figure [AT] is available from the authors, although it should be noted that 
is applicable only for first-order HRC / G800L data which are processed using the same aXe parameters as we adopted 
for our processing, particularly, an extraction window height of ±20 pixels. The methodology is general, however, 
and it would be a straightforward process to re-derive the curve for any combination of processing parameters. Note 
that we used the same methodology (and the same spectrophotometric standards) to derive corrections to the lUE 
Final A rchive calibrat ion iM assa & Fitzpatrick 2000) . In addition, our calibration of the optical and NIR photometry 
(|Fitzpatrick fc Massall200"5b! ). which are used in the modeling of the stellar SEDs, was based on these corrected lUE 
results and, therefore, on the HST white dwarf-based fundamental calibration. With the use of the ACS correction 
derived here, all the data analyzed in this paper are based on an internally-consistent absolute calibration. 

REFERENCES 



Bevington, P. R. 1969, Data reduction and error analysis for the 

physical sciences (New York: McGraw-Hill, 1969) 
Bohlin, R. C. 1996, AJ, 111, 1743 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 
Clayton, G. C, Gordon, K. D., Salama, F., AUamandola, L. J., 

Martin, P. G., Snow, T. P., Whittet, D. C. B., Witt, A. N., & 

Wolff, M. J. 2003, ApJ, 592, 947 
Fitzpatrick, E. L. & Massa, D. 1986, ApJ, 307, 286 
— . 1988, ApJ, 328, 734 
— . 1990, ApJS, 72, 163 
— . 2005a, AJ, 130, 1127 
— . 2005b, AJ, 129, 1642 
— . 2007, ApJ, 663, 320 

Froebrich, D., Murphy, G. C, Smith, M. D., Walsh, J., & Del 

Burgo, C. 2007, MNRAS, 378, 1447 
Gordon, K. 2009, in preparation 

Gosling, A. J., Bandyopadhyay, R. M., & Blundell, K. M. 2009, 
ArXiv e-prints 

Hecht, J., Heifer, H. L., Wolf, J., Pipher, J. L., & Donn, B. 1982, 
ApJ, 263, L39 

^ http:/ /www. stsci.edu/hst/observatory/cdbs/calspec. html 



Johnson, H. L. 1968, Interstellar Extinction (Nebulae and 
interstellar matter. Edited by Barbara M. Middlehurst; Lawrence 
H. AUer. Library of Congress Catalog Card Number 66- 
13879. Published by the University of Chicago Press, Chicago, 
ILL USA, 1968, p.l67). Chapter 5 

Kurucz, R. L. 1991, in NATO ASIC Proc. 341: Stellar Atmospheres 
- Beyond Classical Models, 441 

Lanz, T. & Hubeny, I. 2003, ApJS, 146, 417 

Larsen, S. S., Kiimmel, M., & Walsh, J. R. 2006, in The 2005 HST 
Calibration Workshop: Hubble After the Transition to Two-Gyro 
Mode, ed. A. M. Koekemoer, P. Goudfrooij, & L. L. Dressel, 103 
Larson, K. A. & Whittet, D. C. B. 2005, ApJ, 623, 897 
Markwardt, C. B. 2009, ArXiv e-prints 
Martin, P. G. & Whittet, D. C. B. 1990, ApJ, 357, 113 
Massa, D. & Fitzpatrick, E. L. 2000, ApJS, 126, 517 
Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 
Mathis, J. S. & Wallenhorst, S. G. 1981, ApJ, 244, 483 
Nichols, J. S. & Linsky, J. L. 1996, AJ, 111, 517 



10 



Fitzpatrick & Massa 



Nishiyama, S., Nagata, T., Kusakabe, N., Matsunaga, N., Naoi, T., 
Kato, D., Nagashima, C, Sugitani, K., Tamura, M., Tanabe, T., 
& Sato, S. 2006, ApJ, 638, 839 

Pei, Y. C. 1992, ApJ, 395, 130 

Rieke, G. H. & Lebofsky, M. J. 1985, ApJ, 288, 618 
Spitzer, L. 1978, Physical processes in the interstellaj: medium 
(New York Wiley-Interscience, 1978.) 



Walborn, N. R. 1973, AJ, 78, 1067 

Whittet, D. C. B. 2003, Dust in the galactic environment (Bristol: 

Institute of Physics (TOP) Publishing) 
Zubko, v., Dwek, E., & Arendt, R. G. 2004, ApJS, 152, 211 



IR Extinction 



11 



TABLE 1 

The Target Stars and the new HST/ACS G800L Data 



Star 


Spectral 


ya 




E(B- 




ACS 


Observation 




Type^ 










Datasets 


Dates 


BD+56 517 


B1.5 V 


10. 


.50 


0. 


.51 


J9FH34041, 


J9FH34051 


2005-08-02 


BD+45 973 


B3 V 


8, 


.58 


0, 


.80 


J9FH30041, 


J9FH30051 


2005-09-01 


r>U+44 iUoU 


T3fi TTT 
rSD ill 


y. 


1 o 


U, 


no 


J9FH29041, 


J9FH29051 


zUUu-Uo-zZ 


NGC 1977 #885 


B4 V 


11, 


.33 


0, 


.82 


J9FH15041, 


J9FH15051 


2005-10-10 


HD 46106 


Bl V 


7, 


.92 


0, 


.43 


J9FH38041, 


J9FH38051 


2005-10-19 


HD 292167 


09 HI: 


9, 


.25 


0, 


.72 


J9FH27041, 


J9FH27051 


2005-09-03 


HD 68633 


B5 V 


9, 


.91 


0, 


.49 


J9FH17041, 


J9FH17051 


2006-04-23 


HD 70614 


B6 


9, 


.27 


0, 


.68 


J9FH18041, 


J9FH18051 


2006-09-21 


Trumplcr 14 #6 


Bl V 


11, 


.23 


0, 


.48 


J9FH52041, 


J9FH52051 


2005-11-06 


Trumplcr 14 #27 




11, 


.30 


0, 


.58 


J9FH51041, 


J9FH51051 


2005-08-11 


HD 164740 


07.5 V(n) 


10, 


.30 


0, 


.86 


J9FH08041, 


J9FH08051 


2005-08-19 


HD 229196 


06 ni(n)(f) 


8, 


.52 


1, 


.20 


J9FH23041, 


J9FH23051 


2006-03-14 


HD 204827 


BO V 


7, 


.94 


1, 


.08 


J9FH11041, 


J9FH11051 


2005-08-02 


BD+61 2365 


B0.5 V 


9, 


.19 


0, 


.78 


J9FH42041, 


J9FH42051 


2005-12-22 



^ Spectral types, V magnitudes, and E(B— V) values are taken from Paper V; references for the 
spectral types are given in that paper. 



TABLE 2 

.Ml? Extinction C'riiM'S foh T,\h<:I'T Stars 



k(A-V) 


Star 


K 




H 


J 






AUS87 




(2.18 /im) 




(1.62 ijeca) 


(1.23 ixm) 


(9250 A) 




(8750 A) 


BD-l-56 517 


-2.79 ± 0. 


.09 


-2.65 ±0.12 


-2.34 ± 0.08 


-1.95 ±0. 


.08 


-1.81 ±0.08 


BD+45 973 


-2.47 ± 0. 


.05 


-2.33 ± 0.08 


-2.07 ± 0.04 


-1.58 ±0. 


.05 


-1.46 ±0.05 


BD-l-44 1080 


-2.15 ±0. 


.04 


-2.01 ± 0.06 


-1.80 ±0.04 


-1.38 ±0. 


.04 


-1.26 ±0.04 


NGC 1977 #885 


-4.58 ± 0, 


.08 


-4.24 ±0.10 


-3.47 ±0.06 


-2.43 ±0. 


.06 


-2.21 ±0.05 


HD 46106 


-2.56 ±0, 


.10 


-2.54 ±0.14 


-2.16 ±0.09 


-1.64 ±0, 


.09 


-1.47 ±0.09 


HD 292167 


-2.87 ±0, 


.06 


-2.59 ± 0.09 


-2.31 ±0.05 


-1.78 ± 0, 


.05 


-1.59 ±0.05 


HD 68633 


-3.31 ±0, 


.10 


-3.18 ± 0.14 


-2.76 ±0.08 


-2.17 ±0, 


.08 


-2.00 ±0.08 


HD 70614 


-2.85 ±0, 


.07 


-2.64 ±0.12 


-2.34 ±0.06 


-1.84 ±0, 


.06 


-1.67 ±0.06 


Trumplcr 14 #6 


-4.18 ±0, 


.12 


-3.85 ± 0.14 


-3.27 ±0.10 


-2.32 ± 0, 


.09 


-2.12 ±0.09 


Trumplcr 14 #27 


-3.76 ±0, 


.10 


-3.43 ± 0.12 


-2.87 ±0.08 


-2.05 ± 0, 


,07 


-1.87 ±0.07 


HD 164740 


-4.92 ±0, 


.08 


-4.22 ± 0.09 


-3.49 ±0.06 


-2.48 ± 0, 


,05 


-2.22 ±0.05 


HD 229196 


-2.88 ±0, 


.04 


-2.69 ± 0.05 


-2.34 ±0.03 


-1.74 ±0, 


,03 


-1.59 ±0.03 


HD 204827 


-2.25 ±0, 


.04 


-2.17 ±0.05 


-1.92 ±0.03 


-1.51 ±0, 


.04 


-1.40 ±0.04 


BD-f61 2365 






-2.68 ± 0.07 


-2.30 ± 0.05 


-1.74 ±0. 


.05 


-1.59 ±0.05 



TABLE 2 

NIR Extinction Curves for Target Stars, continued 



k(A-V) 


Star 


ACS82 




ACS77 




ACS72 




ACS67 




ACS 62 






(8250 A) 




(7750 A) 




(7250 A) 




(6750 A) 




(6250 A) 




BD+56 517 


-1.65 + 0. 


.07 


-1.43 ±0 


.07 


-1.18 ±0. 


.07 


-0.92 ±0 


.07 


-0.70 ± 0, 


.07 


BD+45 973 


-1.31 ±0. 


.05 


-1.13 ±0 


.05 


-0.92 ± 0. 


.04 


-0.69 ± 


.04 


-0.44 ± 0, 


.04 


BD+44 1080 


-1.12 ±0. 


.04 


-0.94 ± 


.04 


-0.77 ± 0. 


.04 


-0.60 ± 


.04 


-0.37 ± 0, 


.03 


NGC 1977 #885 


-1.96 ±0. 


.05 


-1.67 ±0, 


.05 


-1.36 ±0. 


.05 


-1.02 ±0, 


.04 


-0.63 ± 0, 


.04 


HD 46106 


-1.31 ±0, 


.09 


-1.16 ±0, 


.08 


-0.96 ±0, 


.08 


-0.76 ± 


.08 


-0.52 ± 0, 


.08 


HD 292167 


-1.43 ±0, 


.05 


-1.26 ±0, 


.05 


-1.05 ±0, 


.05 


-0.82 ±0 


.05 


-0.53 ± 0, 


.05 


HD 68633 


-1.75 ±0, 


.08 


-1.49 ± 


.08 


-1.27±0, 


.07 


-1.06 ± 


.07 


-0.73 ±0, 


.07 


HD 70614 


-1.49 ±0, 


.06 


-1.27 ±0 


.05 


-1.04 ±0, 


.05 


-0.82 ± 


.05 


-0.56 ±0, 


.05 


Trumplcr 14 #6 


-1.88 ±0, 


.08 


-1.62 ±0, 


.08 


-1.36 ±0, 


.08 


-1.05 ± 


.07 


-0.63 ±0, 


.07 


Trumplcr 14 #27 


-1.62 ±0, 


.07 


-1.36 ±0 


.06 


-1.11 ±0, 


.06 


-0.85 ± 0, 


.06 


-0.53 ±0, 


.06 


HD 164740 


-1.95 ±0, 


.05 


-1.67 ±0, 


.05 


-1.36 ±0, 


.04 


-1.04 ±0 


.04 


-0.65 ± 0, 


.04 


HD 229196 


-1.42 ±0, 


.03 


-1.24 ±0 


.03 


-1.04 ±0, 


.03 


-0.82 ± 


.03 


-0.56 ± 0, 


.03 


HD 204827 


-1.26 ±0, 


.04 


-1.10 ±0, 


.03 


-0.92 ±0, 


.03 


-0.71 ±0, 


.03 


-0.43 ± 0, 


.03 


BD+61 2365 


-1.41 ±0, 


.05 


-1.21 ±0, 


.05 


-1.00 ±0, 


.05 


-0.77 ±0, 


.04 


-0.50 ± 0, 


.04 
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TABLE 3 
Results for Power-Law Fits 



Star 




/3-Fixed=' 






I3-Ya.riahle° 








kiR 


/3 R{V) 




kiR 


13 R{V) 












.98 


+0. 
-0 


,06 
.04 


1 


.84 


3. 


03+^ 

.UO_g 


,07 
.07 





.17 


0, 


■°' — 


15 
.10 


2 


.10 


+0. 

-0 


18 

,29 


2, 


95+^ 


13 
.07 





.16 


1.1 




1 


.00 


+0. 
-0 


,04 
.02 


1 


.84 


2, 


73+^' 
■'■^-0, 


,04 
.03 





.07 


0. 


93+^ 


,06 
.06 


2 


.00 


+0. 
-0 


15 

,11 


2, 


68+?' 

.UO_g 


,06 
.06 





.03 


8.8 


BD+44 1080 


0, 


.89 


+0. 

-0, 


,03 
,02 


1, 


.84 


2. 


.391°: 


.03 
.03 


0, 


.29 


0. 


.761°: 


,05 
.04 


2, 


.21 


+0. 

-0, 


15 

,12 


2. 


.281°: 


,04 
.04 


0, 


.02 


63.6 




2 


.13 


+0. 

-0 


,04 
,04 


1 


.84 


4. 


QO + 

■99-0, 


.05 
.07 


2 


.40 


2. 


■97-0, 


.18 
.09 


1 


.26 


+0. 

-0, 


,04 
,08 


5. 


7/1+0. 


,17 
.10 





.64 


14.9 


HD 46106 


1, 


.09 


+0. 

-0, 


07 
,05 


1, 


.84 


2. 


.88t°: 


,09 
.07 


0, 


.23 


1. 


.001°: 


,15 
.10 


2, 


.03 


+0. 

-0, 


,30 
,20 


2. 


.811°: 


,17 
.10 


0, 


.26 


0.4 


HD 292167 


1, 


.18 


+0. 

-0, 


,03 
,04 


1, 


.84 


3. 


.121°: 


,04 
.05 


0, 


.24 


1. 


.281°: 


,10 
.08 


1, 


.66 


+0. 

-0, 


10 

,13 


3. 


.211°: 


,09 
.08 


0, 


.22 


1.4 


HD 68633 


1, 


.35 


+0. 

-0, 


,06 
,04 


1, 


.84 


3. 


.69t°: 


,08 
.07 


0, 


.24 


1, 


.141°' 


,10 
.09 


2, 


.24 


+0. 

-0, 


19 

,16 


3. 


.521°: 


,10 
.10 


0, 


.11 


6.5 


HD 70614 


1, 


.15 


+0. 

-0, 


,04 
,03 


1, 


.84 


3. 


.131°: 


.06 
.05 


0, 


.13 


1, 


.061°' 


,09 
.07 


2, 


.02 


+0. 

-0, 


14 

,21 


3. 


.061°: 


,09 
.07 


0, 


.12 


1.8 


Trumpler 14 #6 


1, 


.90 


+0. 

-0, 


,08 
,05 


1, 


.84 


4. 


.581°: 


.12 
.07 


0, 


.34 


2, 


.341°: 


,20 
.17 


1, 


.44 


+0. 

-0, 


10 

,11 


4. 


.981°: 


,20 
.12 


0, 


.12 


10.7 


Trumpler 14 #27 


1, 


.75 


+0. 

-0, 


,06 
,04 


1, 


.84 


4, 


12+°' 


.08 
.06 


0, 


.35 


2. 


■i6lg: 


,21 
.11 


1, 


.45 


+0. 

-0, 


,09 
,13 


4. 


.481°: 


,20 
.11 


0, 


.03 


53.7 


HD 164740 


2, 


.26 


+0. 

-0, 


,04 
,04 


1, 


.84 


5. 


17+°' 


.06 
.06 


5, 


.51 


4. 


.271°: 


,39 
.24 


0, 


.89 


+0. 

-0, 


,04 
,06 


7. 


.041°: 


,38 
.25 


0, 


.10 


283.1 


HD 229196 


1, 


.23 


+0. 

-0, 


,02 
,02 


1, 


.84 


3. 


.171°: 


.03 
.03 


0, 


.42 


1. 


.291° 


,04 
.04 


1, 


.72 


+0. 

-0, 


,07 
,06 


3. 


.231°: 


,04 
.05 


0, 


.43 


1.0 


HD 204827 


0, 


.86 


+0. 

-0, 


,03 
,02 


1, 


.84 


2. 


.491°: 


.04 
.02 


0, 


.46 


0, 


.731°' 


,04 
.04 


2, 


.24 


+0. 

-0, 


13 

,10 


2. 


.381°: 


,04 
.03 


0, 


.11 


16.5 


BD+61 2365 


1, 


.21 


+0. 

-0 


,05 
,03 


1, 


.84 


3. 


.141°: 


,05 
.05 


0, 


.09 


1, 


.381S' 


,34 
.17 


1, 


.61 


+0. 

-0, 


,28 
,34 


3. 


.301°: 


,32 
.16 


0, 


.05 


3.7 



^ "Universal'' value of /3 — 1.84 assumed in these fits. Reduced-^ values are shown, corresponding to 5 degrees 
of freedom. Fits for the sight line towards BD+61 2365, for which no i<'-band data are available, have 4 degrees of 
freedom.^ Power-law exponent /3 taken as a free parameter in these fits. Reduced-x^ values are shown, corresponding 
to 4 degrees of freedom. Fits for the sight line towards BD+61 2365 have 3 degrees of freedom.^ The F-statistic, 
which tests the need to include ^ as a free parameter in the power law fits jBevingtonll 19691) . 



TABLE 4 

Results for Fits Using Equation {sjl^ 



Star 


a 


R{V) 


X^ X^ 








(/3-Variablo) 



BD+56 517 


2, 


OQ + 0.12 

■«9-0.15 


2. 


.86: 


BD+45 973 


2, 


49+0.08 
■^^-0.10 


2. 


.64: 


BD+44 1080 


2, 


47+0.06 
■^'-0.10 


2. 


so: 


NGC 1977 #885 


1, 


7O+0.05 
■ '''-0.06 


5. 


.54: 


HD 46106 


2, 


4^+0.20 
■^''-0.25 


2. 


.77: 


HD 292167 


2, 


40+0-10 


3. 


.04: 


HD 68633 


2, 


=4+0.09 
■0^-0.15 


3. 


.51; 


HD 70614 


2, 


48+0-09 
■48-0.10 


3. 


.02; 


Trumpler 14 #6 


1, 


94+0-08 
■''^-O.ll 


4. 


.80: 


Trumpler 14 #27 


1, 


fil +0.08 


4. 


.45: 


HD 164740 


1, 


61+°-°'* 
■"-'--0.05 


6. 


.oo: 


HD 229196 


2, 


0^+0.05 
■■^^-0.05 


3. 


.09: 


HD 204827 


2, 


7t:+0.07 
■ '''-0.09 


2. 


.35: 


BD+61 2365 


2, 


OO+O.ll 

■^"-0.12 


3. 


.ii; 



+0.07 
-0.07 
+0.06 
-0.04 
+0.05 
-0.03 
+0.12 
-0.11 
+0.12 
-0.08 
+0.06 
-0.07 
+0.11 
-0.08 
+0.07 
-0.06 
+0.21 
-0.16 
+0.19 
-0.11 
+0.15 
-0.12 
+0.04 
-0.04 
+0.04 
-0.02 
+0.11 
-0.09 



0.20 


0.14 


0.13 


0.23 


0.30 


0.69 


0.38 


0.87 


0.17 


0.33 


0.21 


0.22 


0.38 


0.60 


0.22 


0.48 


0.10 


0.24 


0.18 


0.43 


0.15 


0.20 


0.75 


1.27 


0.13 


0.24 


0.08 


0.19 



^ The final column is foi' fits using the /3- Variable model applied 
to the same wavelength range used in the fits employing equation | f5l . 
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Fig. 1. — Final processed ASC/HRC G800L spectra for the 14 program stars. A Kurucz ATLAS9 model with Teff= 20000 K is shown 
near the top of the figure for comparison. Weak Ha and H I Paschen absorption lines may be present in our spectra, although they would 
be close to the detection limit. All the prominent "bumps and wiggles" seen in the ACS spectra are instrumental in origin. 
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Fig. 2. — Multiwavelength SEDs for the program stars. In the UV (A < 3000 A), low-resolution lUE spectrophotometry is shown as 
small circles. In the optical (3000 A < A < 6000 A), Johnson UBV, Stromgren uvby, and Geneva UB^B-z V\G photometry are indicated 
by the circles, triangles, and diamonds, respectively. In the NIR, the small circles are the new ACS observations and the large circles are 
the 2MASS JHK photometry. The SEDs have been offset by arbitrary amounts for display purposes. The dashed lines are not model fits 
to the SEDs, but were added to guide the eye and to help distinguish the datasets for each star. 
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Fig. 3. — NIR extinction curves for six of the 14 program sight Unes. The different symbols show the fc(A — y) for the 2MASS JHK 
photometry (large circles), ACS spectrophotometry (open and filled squares), and Johnson V magnitudes (triangles), as listed in Table 2. 
The ACS dat a we re binned into a set of seven magnitudes before computing the extinction. The unbinned ASC data are shown as small 
circles (see S I3.1|l . The dashed curves are the fits to the JHK and ACS data using Equation l(2} with /3 = 1.84 ("/3-Fixed" results in 
Table 3). The solid curves are fits using Equation J21I with /3 as a free parameter ( "/3— Variable" results in Table 3). For both sets of fits, 
only ACS data longward of 7500 A (A"'^ < 1.33 were included (filled squares). Shorter wavelength data (open squares) are not 

consistent with a single power-law representation. The dotted curves are /3-Variable fits to the ACS data only. The curves are organized 
in order of increasing Riy). 
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Fig. 3 (coNT.). — 



Continuation of Figure |3] for six more sight lines. 
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Fig. 4. — The power-law exponent f3 versus E{B — V) (top panel) and R{V) (bottom panel) for the "/3— Variable" fits from Table 3. The 
dotted horizo ntal lines show the value /3 = 1.84. The filled triangle in the lower panel represents the sight line to the star p Oph from 
IMartin fc Whittetiil99a The I-ct error bars are based on fits to the Monte Carlo simulations of the data and the noise model. 
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Fig. 5. — Observed NIR extinction curves compared to fits determined from Equation JSjl. The data and symbols are tlie same as in 
Figure [3] All the ACS data were included in the fits (thus all the large squares are shown as filled). The dotted portions of the curves 
indicate extrapolations of the fitting function beyond the range of the data. The sight lines are organized in order of increasing R{V). 




Fig. 5 (coNT.). — 



Continuation of Figure [5] for six more sight lines. 
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Fig. 5 (CONT.). — Continuation of Figure|5]for the two highest R{V) sight hnes. 
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Fig. 6. — Plot of the a and -R(V) from Table 4. The l-cr error bars are based on fits to the Monte Carlo simulations of the data and the 
noise model. 
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Fig. 7. — Comparison of extinction curve parameters derived from observed and theoretical extinction curves, cimwsi is the exponent 
that results from a simple power-law fit (see eq. [2]) to the IR extinction at = 0.00, 0.29, and 1.11 fira~^. It provides a measure of 
the shape of the curves for A > 9000 A. The crosses show ctMWSl ^'^'^ calculated from the lMathis &: Wallenhorsti }1981') theoretical 

extinction curves for a variety of size distributions and abundance ratios (see text). The small, medium, and large crosses correspond to 
(C/Si)tj„st = 3.25, 6.5, and 13.0, respectively. The filled circles show the same quantities derived from the best-fitting NIR curves shown 
in Figure [5] This figure shows that some combinations of model dust abundance and size distribution are not favored by nature. Most 
notable, are the crosses in the lower left of the figure, which does not contain observed points. The models in this region correspond to 
size distributions with small upper cutoffs for silicates {a^ = 0.25 ^m) and large upper cutoffs for graphite {a+ = 0.40 or 0.50 fim). The 
IMathis fc WallenhorsH II1981I ) models are not useful for addressing the shapes of the smaller R{V) sight lines. 
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Fig. 8.— Plot of the observed k{K - V) = E{K - V)/E(B - V) versus the R{V) derived from fits to Eq. The 1-a error bars 

illustrate the correlated errors in the two quantities. The dashed line gives the often-used relation R(V) = —1.1 i which is based 

on van de Hulst's theoretical extinction curve No. 15 (see, e.g.. IJohnson|[T968l ). An unweighted least squares linear fit to the data is shown 
as the solid line, R{V) = -1.36 %\%Zv) " 
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Fig. A2. — Verification of the ACS/HRC/G800L calibration correction from Figure fAll The main figure shows the observed SED 
of HD 47839, including low-resolution lUE spectrophotometry (small circles) in the UV (A < 3000 A), UBV (large circles), Stromgren 
uvby (triangles), and Geneva UB1B2V1G (diamonds) photometry in the optical (3000 A < A < 6000 A), and RIJHK photometry in 
the NIR (A > 6000 A). The smooth curve is a 37,000 K TLUSTY model atmosphere IILanz fc Hubenvl [2053) fit to the data from (e.g., 
IFitzpatrick fc Massa 2005a b). The inset illustrates the effect of the ACS/HRC/G800L correction curve shown in Figure ED It shows the 
best fit model from the main figure (continuous curve) with the corrected and uncorrected G800L data (points) overplotted. For A < 9500 A 
the re-calibrated data mimic the shape of the model SED, but lie 1% above it, consistent with the scatter among the spectrophotometric 
standards. For A > 9500 A, the fiuxes are clearly much less reliable 



